Chirality sensitive effect on surface states in chiral p-wave superconductors 
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We study the local density of states at the surface of a chiral p-wave superconductor in the presence 
of a weak magnetic field. As a result, the formation of low-energy Andreev bound states is either 
suppressed or enhanced by an applied magnetic field, depending on its orientation with respect to 
the chirality of the p-wave superconductor. Similarly, an Abrikosov vortex, which is situated not too 
far from the surface, leads to a zero-energy peak of the density of states, if its chirality is the same as 
that of the superconductor, and to a gap structure for the opposite case. We explain the underlying 
principle of this effect and propose a chirality sensitive test on unconventional superconductors. 
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Much attention has been paid to unconventional su- 
perconductors, because they can exhibit a sign or gen- 
eral phase change of their gap function as a function 
of momentum. This property induces many intriguing 
phenomena, which can be observed directly by so-called 
phase sensitive experiments providing powerful tools to 
test the symmetry of the gap function 1] . One important 
consequence of the sign change of the gap function is the 
possible existence of Andreev bound states at the surface 
of the superconductor [1,[1,|3|- The formation of Andreev 
bound states increases the local zero-energy quasiparti- 
cle density of states (DOS) at the surface, leading to a 
pronounced zero-bias conductance peak in the tunneling 
conductance observable both in singlet d-wave supercon- 
ductors like the cuprates and in tripletp-wave supercon- 
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ductors such as Sr2Ru04 
For the case of d-wave superconductors it is well-known, 
that an applied magnetic field or an applied electric 
current result in a split of this zero-bias conductance 
peak, since the zero-energy spectral weight of the bound 
states is effectively Doppler shifted towards higher ener- 
gies d, Hi, 15, 16 1 . The same effect also appears for an 
Abrikosov vortex, which is pinned not too far from the 
boundary. Here, the zero-energy DOS is suppressed in a 
shadow- like region 'behind' the vortex (itI. [lit. 

Regarding the chiral p-wave superconducting phase as 
it is likely realized in Sr2Ru04, a further aspect ap- 
pears. The chiral p-wave state characterized by the vec- 
tor d(k) = (0, 0, kx iiky) breaks time reversal symmetry 
[3, S, EH, [2I 0|. In this Letter we will show that 
the influence of an external magnetic field on the sur- 
face density of states is selective for the chirality. The 
quasiparticle density of states at the surface increases 
or decreases depending on the relative orientation of the 
applied magnetic field and the chirality. Similarly, we 
find that the influence of a vortex on the surface states 
depends on the orientation of vorticity with respect to 
chiralty. These characteristic effects could open an al- 
ternative way to chirality sensitive probes, in contrast to 



phase or spin sensitive setups, which have been inten- 
sively used in the scientific community already. 

For our calculations, we use quasiclassical Eilenberger 



theory of superconductivity [2J, |25| in the so-called 
Riccati-parametrization ^26'| , which allows to achieve nu- 
merically stable solutions for the quasiclassical propaga- 
tors (and thus also for the DOS) in spatially nonhomo- 
geneous systems. Concretely, in our case we consider a 
superconducting half space a; > exhibiting a gap func- 
tion of chiral p-wave symmetry. The surface at a; = 
is included in our calculations in a straightforward way 
by specific boundary conditions [13, [28, 29]. We assume 
a cylindrical Fermi surface for the superconductor with 
the symmetry axis pointing along the z-direction, so that 
the Fermi velocity can be parametrized by the polar an- 
gle 9 via vp = vp {x cos -\- y sin 6) . For a given point 
Tq in space and angle 9 parametrizing the Fermi surface, 
a quasiclassical trajectory is then defined according to 
r{x') — ro -\- x' Vp. Along such a trajectory the Eilen- 
berger equations can be transformed into 2x2 matrix 
differential equations in spin space, which are of the Ric- 
cati type and can be solved much easier [29j . In our case, 
we deal with a one-component gap-function, so that the 
corresponding Riccati equations take the simpler form 

Tivpdx'a{x') -I- [2e„ 4- A''a{x')] a(x') - A = 
nvpdx,b{x')-[2in-^Ab{x')]b{x')-hA^ = (1) 

for the two scalar coherence functions a and b. Here, 
iin — ien-^-vp ■ - A denotes Matsubara frequencies which 
are shifted due to the presence of a magnetic vector po- 
tential A. The pairing potential A can be factorized in 
the following form 



A{r,e) = Aoexp(i0)^'(r). 



(2) 



Here, ^' denotes a factor which covers the spatial depen- 
dence of the pairing potential in general. Since we are 
only interested in the main qualitative aspects of the lo- 
cal DOS, namely, if the zero-energy spectral weight at 
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the surface is suppressed or increased, we may take the 
modulus of \1/ to be constant 0, 
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For the cal- 
culation of physical properties, the Riccati equations ([J) 
have to be integrated numerically using proper starting 
values in the bulk. The local DOS, which is already nor- 
malized to the DOS in the normal state, is then achieved 
by an integration over the Fermi surface. In terms of the 
coherence functions a and b, we have 



N{ro,E) = 
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where E denotes the quasiparticle energy with respect to 
the Fermi level and 5 is an effective scattering parameter 
that corresponds to an inverse mean free path. For all 
numerical calculations, we fix this value as 5 — O.IAq. 

In order to study the basic effect of chirality we con- 
sider a magnetic field applied along the z-axis at the sur- 
face, represented by a nearly homogeneous vector po- 
tential A. We choose the real gauge, i.e. the spatially 
dependent part ^' of the pairing potential is taken to be 
real. Since we additionally assumed a spatially constant 
modulus of the pairing potential, it is possible to get an- 
alytical solutions for the coherence functions a and b in 
this case, which also allows to examine the corresponding 
behaviour of the local DOS analytically. Directly at the 
surface, we get 



N{E) = 2Re 
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(4) 



with ttin = sAoe^^^~ ' , bout — sApe" and the abbrevia- 
tion s = l/{in + i/e^ -I- Aq). Furthermore, (...) denotes 
angular averaging, which we may restrict to outgoing an- 
gles —tt/2 < 9 < tt/2 only. This directly yields 



N{E) = 2Re 
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Expanding the zero-energy DOS in orders of the vector 
potential A, we obtain in the clean limit of 5 — *■ 0^ 



N{E = 0) = 1 + -^Ay 



(6) 



Physically, this result displays the influence of a Doppler 
shift due to a superfluid velocity on the local quasipar- 
ticle spectrum. In terms of chiral surface states 30,, ..31] , 
the Doppler shift leads to a change of the slope of the 
quasiparticle dispersion [(.(ky) = Aoky/kp], which has a 
direct impact on the corresponding DOS. As is seen from 
the result on the right-hand side of Eq. ([6]), the term of 
the vector potential, which belongs to the direction per- 
pendicular to the chirality of the p-wave superconductor, 
survives in linear order. This is in contrast to other super- 
conductors like S-, d- or p-wave superconductors without 
chirality, where we obtain similarly 



N{E = 0) = C+ {F{d) sin 0) Ay + 



(7) 
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FIG. 1: (color online) N{E = 0) in a chiral p-wave super- 
conductor in the presence of a vortex, which is situated at a 
distance of xv = 2^ from the surface. The vortex and the 
p-wave state have (a) the same chirality and (b) the oppo- 
site chirality. Clearly, there is a remarkable difference in the 
quasiparticle spectral weight at the surface {x = 0) behind the 
vortex, showing a strong increase in (a) and a suppression in 
(b). 



with a constant C and a function F which satisfies 
F{9) = F{—9). Thus, after angular averaging, the coeffi- 
cient of the linear term vanishes in these cases, reflecting 
the presence of inversion symmetry with respect to the 
x-y plane. Since an applied magnetic field is related to 
the vector potential as Bz = dxAy, the zero-energy DOS 
in the chiral p-wave superconductor depends on the sign 
of the magnetic field. Applying a weak magnetic field 
along the chirality direction suppresses the zero-energy 
bound states, while applying it in the opposite direction 
leads to a zero-energy peak of the surface DOS. It is im- 
portant to realize that the derived Eq. ^ qualitatively 
implies this chirality sensitive effect also for the case of 
a more general vector potential. Especially, this chirality 
effect is remarkable in the presence of a vortex near the 
surface. 

As a next example of the chirality effect, we study the 
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FIG. 2: (color online) Local DOS at the point x = y — 
for different vortex to boundary distances xv as a function of 
energy. Again, p-wave state and Abrikosov vortex have the 
same chirality in (a) and the opposite in (b), resulting in peak 
and dip structures, respectively, around zero energy. 



case of a single Abrikosov vortex line parallel to the z- 
axis. For such a vortex in the bulk of a chiral p-wave 
superconductor, the interplay between vorticity and chi- 
rality has been intensively studied already |32l |. In the 
present work, however, we focus on important surface ef- 
fects, which appear at the boundary due to the presence 
of the vortex. We assume the vortex to be at a distance 
Xv from the boundary. With y denoting the coordinate 
along the boundary, the vortex shall sit at yv' = 0. For 
convenience, we set the vector potential A = and in- 
corporate the vortex by a properly chosen phase factor 
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tion z = X + iy) is given by 17|, [l 



instead, which (in standard complex nota- 



i^{r) _ Z- Zy / Z- Zy 
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Here, the first factor is the phase of a single vortex at po- 
sition zy, while the second factor corresponds to a virtual 
antivortex placed at the mirrored position zy — —xy, 
which ensures the correct implementation of boundary 
conditions. In the following, we consider two different 
cases: The chiralities of the p-wave state and the vortex 
are the same (a), or opposite (b). The latter case is im- 
plemented by a complex conjugation of the phase factor 



given in Eq. ([8]), replacing the vortex by an antivortex 
and vice versa. Our results for the corresponding local 
zero-energy DOS near the surface of the p-wave supercon- 
ductor are shown in Fig. [TJ The Abrikosov vortex posi- 
tion is fixed at a distance of xy = 2^ from the boundary 
with ^ — hvp/^o denoting the coherence length. Apart 
from zero-energy bound states in the vortex core, there 
are also bound states at the surface of the p-wave su- 
perconductor. Far away from the vortex, these surface 
bound states have the spectral weight of N{E = 0) = 1 
[cf. Eq. dS])]. However, as can be seen quite clearly in 
Fig. [H the local DOS drastically changes in a shadow- 
like region behind the vortex. If vortex and p-wave state 
have the same chirality, the bound states are strongly 
enhanced (a), for opposite chirality they are suppressed 
(b) . The latter suppression resembles a similar effect ap- 
pearing in d-wave superconductors |17|, |18|. 

In Fig. [21 we show the local DOS at the point x — 
y = for different vortex to boundary distances xy as a 
function of energy. Around zero energy, we find a sharp 
peak or dip structure, respectively, again depending on 
the chirality. These structures get less pronounced, when 
the vortex distance is increased, nevertheless they per- 
sist. Moreover, the quasiparticle spectrum starts to ex- 
hibit some kind of mirror symmetry around the value 
N{E) = 1 for the two different chiralities. Note that 
this is in qualitative agreement with Eq. ([6]) since after 
a transformation to the real gauge we eventually have 
Ay > due to the vortex, leading to the strong increase 
of Andreev bound states at the surface for the same chi- 
rality. The result for opposite chirality is obtained due 
to Ay < for the antivortex, accordingly. It is worth 
noting that the modification of the surface quasiparti- 
cle states due to the presence of vortices has an effect 
on the force acting on vortices near the surface. An in- 
crease of the DOS leads to a repulsion of the vortex from 
the boundary towards the bulk, whereas a decrease re- 
sults in an attraction towards the boundary. Thus, in 
both cases the Bean-Livingston barrier would be modi- 
fied, which influences the escape and entrance of vortices 
to the superconductor [s^, . 

Our results allow us to propose a chirality sensitive 
test on superconductors, based on well-estabhshed ex- 
perimental techniques, which are capable of indicating 
the weight of Andreev bound states, for example tunnel- 
ing conductance experiments or low-temperature scan- 
ning tunneling spectroscopy (35^ . For a chiral supercon- 
ductor, it is expected to observe a suppression of the 
zero-energy DOS at the surface, when a weak magnetic 
field is applied parallel to the chirality. Inverting the 
field, however, leads to an enhancement of the DOS. In 
this way chirality could be detected. This unusual rever- 
sal effect does not appear in non-chiral superconductors. 
Moreover, the experiment allows us to detect the chiral- 
ity and possibly even the domains of different chirality, 
if domain walls reach the surface. 
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In summary, we have studied the DOS at the surface 
of a chiral p-wave superconductor in the presence of a 
magnetic vector potential due to, for example, an ap- 
plied magnetic field or an Abrikosov vortex. We clarified 
that the weight of low-energy surface bound states gets 
either suppressed or increased, depending on the orien- 
tation of both chirality and magnetic field. Due to this 
chirality-sensitive effect on the Andreev bound states, a 
setup to test the chirality of an unconventional supercon- 
ductor could be accessible experimentally. Furthermore, 
for vortices this effect also has a chirality selective influ- 
ence on the Bean-Livingston barrier which could give rise 
to a different escape rate of vortices from the two kinds 
of chiral domains. 
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